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Abstract Measurements of the scalar dissipation rate 
are performed in an impulsively started gas jet, using 
planar laser induced fluorescence. The measurements 
are well resolved spatially. The deteriorating effect of 
experimental noise on this experiment is treated with 
a Wiener filter, which is shown to be applicable to 
this large-scale inhomogeneous flow. The accuracy of 
the scalar dissipation rate is within 20%, as determined 
from an explicit calculation of the filtering errors. The 
residual fields that remain after the filtering are anal¬ 
ysed in detail and their statistical properties show that 
these resemble white noise to a good approximation. 
The level of corrections is minimal for the scalar field 
but it is of the order of 40% for the scalar dissipation 
rate. An examination of the filtering operation using 
modeled spectra and the measured spatial resolution 
shows that the Wiener filter produces errors in the es¬ 
timate of the scalar dissipation rate ^ 30%, for Taylor- 
scale Reynolds number up to 1000. The implications of 
this modelling are discussed with respect to common 
experimental situations and point out the relative mer¬ 
its of improving the spatial resolution as compared to 
improvements in the signal to noise ratio. 
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1 Introduction 


The dissipation rate of scalar fluctuations in a turbulent 
flow, X, is defined as 


x = D 


{dQ/dxf' + {dQ/dyf' + {dQ/dzf 


( 1 ) 


where C is a scalar and D is its diffusivity. It is of fun¬ 
damental interest in mixing between fluids, [Sr eenivasan 
(|1991|) ; [Warhaftj ([2QQQ|) ; [Dimotakisj (|2QQ5|) , and is a mea¬ 


sure of the degree of active mixing that is in progress in 
a turbulent flow. For example, in homogeneous turbu¬ 
lence, the scalar fluctuations are ’smoothed out’ by the 


scalar dissipation rate, Tennekes and Lumley (1972), 


and as a consequence its increase leads to faster molec¬ 
ular mixing between fluids. Considerable attention has 
been paid to the measurement of the scalar dissipation 
in variety of flows, which has served to establish basic 
and common characteristics of this quantity in a va¬ 


riety of flows. For example, Antonia and Mi (1993a); 


Anselmet et al (1994) used cold wires to measure the 


dissipation of temperature fluctuations (considered as a 
passive scalar) in axially symmetric jets an d [Anselmet 
and Antonia ( [1985 ); Antonia and Browne (1983) used 
the same technique in plane jets. In [Kailasnath et al 
(1993), a variety of flows and techniques were used to 


identify common characteristics of the scalar dissipa¬ 
tion rate in both gaseous and liquid flows. Furthermore, 


Dowling (1991) and Pitts et al (1999) used point and 


line, respectively, Rayleigh scattering to measure time 
series of the concentration and its dissipation in round 
gas jets. The spatial structure of the scalar dissipation 
field was demonstrated, in liquid and gaseous round 
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jets, by Buch and Dahm (1996) and Buch and Dahm 


(1998), respectively. They used planar laser induced flu¬ 
orescence and planar Rayleigh scattering, respectively, 
to record 2-diniensional images of instantaneous dissi¬ 
pation structures; measurements of these structures in a 
gaseous planar jet were performed by |Su and Clemens 


(2003). Additionally, 2-dimensional Raman scattering 


was used, Namazian et al (1988); Schefer et al (1994) 


to measure the scalar dissipation rate in an axisym- 
metric gaseous jet in order to elucidate the role of large 
scales on the jet development. Three-dimensional mea¬ 
surements of the scalar dissipation rate have, also, been 
attempted in gaseous and liquid jets and wakes, for ex¬ 


ample by Prasad and Sreenivasan (1990). Other flows 


have also received attention, for example grid turbu¬ 


lence, Tong and Warhaft (1994), where passive temper¬ 


ature fluctuations and their dissipation were measured 
using cold wires. Finally, the scalar dissipation rate in 
opposed jet flows was measured using either pairs of 


cold wires, Sardi et al (1998), or planar laser induced 


fluorescence of acetone, 


Krawczynski et al (2006a) 


Furthermore, in the case of non-premixed flames, 
chemical reaction happens only after molecular mixing 
between fuel and oxidiser, so the scalar dissipation rate 
is of prime importance in combustion as well, where 
it arises naturally in the expression of the mean reac¬ 
tion rate, Bilger (1976). Additionally, in the flamelet 


modelling approach of non-premixed turbulent flames 
Williams (1975); Peters ( 1984| ) the scalar dissipation 
rate of the mixture fraction (a normalised fuel concen¬ 
tration) provides the link between chemistry and flow 
field. There is evidence as well, e.g. [Mastorakos et al 


(1997), that the scalar dissipation rate is a key quan¬ 


tity controlling the autoignition of two, initially sep¬ 
arate, gases. Despite their considerable difficulty, mea¬ 
surements of the scalar dissipation rate of mixture frac¬ 
tion or temperature in non-premixed flames have been 
performed. For example, Rayleigh scattering has been 
performed to acquire point or 2-dimensional images of 
temperature and its dissipation in turbulent jet flames. 


e.g. Everest et al ( 

1995); Effelsberg and Peters (1988); 

Kaiser and Frank 

to 

o 

o 

), and it was, also, combined 


with planar laser induced fluorescence to measure fuel 


concentration, Sutton and Driscoll (2002). Many mea¬ 


surements of the dissipation rate of mixture fraction in 
turbulent jet diffusion flames have been performed us¬ 
ing line Raman scattering, combined with laser induced 


and Mansour 


fluorescence; for example, Nandula et al (1994); Chen 


The scalar dissipation rate is a small scale quantity, 
so the measurement probe volume should ideally be of 
the order of the Kolmogorov or Batchelor length scales, 
George and Hussein (1991). However, at these small 


resolutions experimental noise can become comparable 
with the measured signal. Many studies have shown the 
effect of experimental noise and spatial resolution on 
the measurement of the scalar dissipation rate using 
either analytical techniques, computational results or 


experimental data, for example Wyngaard (1971); An- 

tonia and Mi 

(1993b); 

Mi and Nathan (2003 

); Burattini 

et al (2007); 

Ghandhi 

(2005). It is generally true that 


in ideal, noise-free measurements, the magnitude of the 
scalar dissipation increases with increasing spatial reso¬ 
lution and converges to the true value for well-resolved 
measurements. Similarly, the experimental noise tends 
to increase the magnitude of the scalar dissipation with 
increasing spatial resolution. However, it is known that 
as the measurement probe volume becomes very small, 
this increase is unbounded and faster than the increase 
due to resolution effects; it is possible to measure many 
times the actual average value of the dissipation due 
to noise effects. So, under-resolved, but noisy, measure¬ 
ments can give the ’correct’ value of (y). Therefore, 
measurements of the mean scalar dissipation rate have 
to be checked, and possibly corrected, for noise and 
spatial resolution effects and measurements of the in¬ 
stantaneous scalar dissipation require post-processing 
of the scalar measurement, typically in the form of dig¬ 
ital filtering for noise removal. 

It is possible to check that the mean scalar dis¬ 
sipation rate has been measured correctly by either 
using theoretical arguments (mentioned below) or by 
verifying the balance of the scalar variance transport 


equation, e.g. Antonia and Mi (1993a); Bilger (2004). 


Furthermore, resolution correction factors for the mean 
scalar dissipation, in homogeneous and isotropic turbu¬ 
lence, rely in averaging a three-dimensional spectrum 
of the scalar fluctuations, either using a modelled form 


Wyngaard (1968 1971), or using measured one-dimensional 


spectra Antonia and Mi| (1993b). These results show 
that the spatial resolution required to measure the mean 
scalar dissipation rate within 10% of the true value is 
2 ^ 3r]B {vb being the Batchelor length scale), when 
cold wires are used for the measurement of passive tem¬ 
perature fluctuations. 

Various methods have been devised to calculate the 
mean scalar dissipation rate from measurements of the 
underlying scalar. In a slightly heated axially symmet- 


(1997) measured hydrogen jet flames, whereas ric jet Antonia and Mi (1993a) employed a pair of cold 


Starner et al (1994); Karpetis and Barlow (2002); Geyer 


et al (2005) measured the dissipation rate in hydrocar¬ 


bon flames. 


wires and used two techniques for the determination of 
the mean scalar dissipation rate. The measured temperature- 
derivative spectra were used to devise correction fac- 
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tors for the mean scalar dissipation rate, by adapting 


the above mentioned analysis Wyngaard (1971); fur¬ 
thermore, a Taylor series expansion of the cross corre¬ 
lation coefficient for small distances was extrapolated 
to zero separation, giving an estimate of the mean dis¬ 
sipation rate. The two techniques agreed very well and 
the accuracy of the mean scalar dissipation rate was 
determined by showing that it balances the scalar vari¬ 
ance transport equation. The cross correlation coeffi¬ 
cient technique can also be used for non-homogeneous 
flows by including higher order terms in the Taylor se¬ 


ries expansion, Krishnamoorthy and Antonia (1987). A 


method that relies on the second order structure func¬ 
tion estimates the noise variance by plotting the mean 
squared scalar difference for various separations. An ex¬ 
trapolation of this plot to zero separation gives a mea¬ 
sure of the noise variance, as was done in a heated round 


jet Anselmet et al (1994) and in a heated opposed jet 


flow Sardi (1997), using pairs of cold-wires. Then, the 
second order structure function, minus the noise vari¬ 
ance, is divided by the squared separation and plotted 
as a function of the separation. A final extrapolation to 
zero separation gives the mean scalar dissipation rate. 
Although the above methods measure accurately the 
mean scalar dissipation rate, it is always desirable to 
ensure a high signal-to-noise ratio. It should be noted 
that these methods are mainly applicable to homoge¬ 
neous and isotropic turbulence. 

Other techniques, for the measurement of the mean 
scalar dissipation rate, rely on over-sampling or redun¬ 
dant scalar measurements. For example, [Kaiser and Frank] 
(2007) measured the temperature signal in a turbulent 


non-premixed flame using 2-dimensional laser Rayleigh 
scattering. The fact that noise is uncorrelated between 
adjacent rows (or columns) of the image data allows 
the calculation of power spectra by interlacing between 
consecutive rows. The noise level was reduced in the 
scalar derivative spectra, from which the mean scalar 
dissipation can be calculated. Also, in experiments in 
a turbulent jet flame using single-point laser Rayleigh 


scattering Wang et al (2005) used either over-sampling 


or two different probe volumes to measure the tempera¬ 
ture fluctuations at the same spatial location. The noise 
variance was then identified and subtracted from the 
measurement of the mean scalar dissipation and cor¬ 
rections of up to 90% of the uncorrected dissipation 
were obtained. Finally, the kinetic energy dissipation 
rate was calculated at two different resolutions from a 


single set of data Tanaka and Eaton (2007) and, by com¬ 


bining the two measures, the correct mean dissipation 
rate was inferred (this technique can be readily adapted 
to the measurement of the scalar dissipation rate). Since 
these techniques do not rely on the assumptions of ho¬ 


mogeneity and isotropy, they seem broader in scope for 
the determination of the mean scalar dissipation rate. 
However, they do not provide accuracy estimates for the 
resulting measured mean dissipation rate (the method 
in Tanaka and Eaton (2007) does provide a measure of 


its accuracy). 

The measurement of the instantaneous scalar dissi¬ 
pation rate is more restrictive in terms of signal-to-noise 
ratio, while being equally important. So, it generally re¬ 
quires the application of digital filtering, or other post¬ 
processing, to the raw scalar measurements. In measur¬ 
ing the passive temperature fluctuations in a jet flow 
using cold wires. Mi et al (2005) used an iterative pro¬ 
cedure to filter the raw data, by continuously adjust¬ 
ing the cut-off frequency of a digital filter. The algo¬ 
rithm used Kolmogorov scalings and converged to the 
digital cut-off frequency when this value and the up¬ 
dated Kolmogorov frequency were the same within a 
very small bound. In jet flame experiments Kaiser and] 


Erank (2007) obtained the instantaneous scalar dissipa¬ 


tion rate by filtering the instantaneous 2D temperature 
images using gaussian kernels, whose widths were dif¬ 
ferent at each of the two spatial directions and, also, 
varied with the position in the flow. Recently, jCai and 


Tong (2009) proposed a method that relied, also, on 


the measured scalar field and the expectation that the 
so-called sub-grid scale variance of the mixture fraction 
is small at well-resolved, low-noise portions of the data. 
In this way, the noise variance could be obtained and 
used to correct under-resolved data. In other experi- 
ments, [Miller and Dimotakis (1996); Krawczynski et al 


(2006a) used an optimal Wiener filter approach, that is 
also adopted in the present experiment, to filter Laser 
Induced Eluorescence (LIE) measurements in a round 
jet and in a partially stirred reactor, respectively. Ei- 
nally, Markides and Mastorakos (2006) used a filtering 
method that relies on a wavelet decomposition of the 
raw scalar measurements. 


The above exposition shows that in the case of ho¬ 
mogeneous and isotropic turbulence, the signal process¬ 
ing involved in measuring the mean scalar dissipation 
rate can be checked against theoretical considerations 
to ascertain the accuracy of the measurement. In a more 
general flow configuration, concerning also the measure¬ 
ment of the instantaneous scalar dissipation rate, the 
signal processing takes a central role in determining 
the scalar dissipation rate, so it is important to fully 
analyse the post processing procedure. Eor example, 
filters usually rest on certain assumptions about the 
underlying signal and it is indeed possible to check the 
consistency of the filtering operation. This is necessary 
because the effects of filtering, for example the level of 
corrections, the residual fields and their characteristics. 
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are not usually demonstrated for any given filtering op¬ 
eration. Furthermore, many of the existing methods for 
correcting the mean scalar dissipation rate do not pro¬ 
vide corrections for the instantaneous scalar fields. 

In the present paper we develop an approach to 
address the issues above by using the Wiener filtering 


theory Miller and Dimotakis (1996). We show how to 


explicitly calculate the errors that arise in the filter¬ 
ing process and so give an accuracy estimate for the 
mean scalar dissipation rate. We, also, present in de¬ 
tail the effects of filtering on the measured scalar fields. 
It is important to note that the Wiener filter provides 
both a noise and a resolution correction (since it makes 
use of the Point Spread Function of the measuring sys¬ 
tem) for the instantaneous scalar fields. We show that 
the Wiener filter is rather insensitive to the modelling 
involved, while at the same time is easily applied; it 
is, however, a ’global’ method, since it relies on scalar 
power spectra. In order to verify that the noise and 
the resolution corrections are accurate, we compare the 
mean scalar dissipation obtained from the Wiener filter 
with results from another technique [Tanaka and Eaton] 
(2007), which has been shown to correct for noise and 


resolution effects as well as being second order accu¬ 
rate in the spatial resolution (thus giving an estimate 
of the accuracy). However, we choose the Wiener fil¬ 
ter because we can filter instantaneous images and ob¬ 
tain the distributions of the scalar dissipation rate and 
its unconditional and conditional (on the scalar value) 
statistics. 

We use a transient, incompressible jet flow in our 
measurements, since this flow is of technological im¬ 
portance, e.g. gaseous fuel injection in engines, and, 
also, has received limited attention in terms of mix¬ 
ing measurements. It is more common to report veloc¬ 
ity measurements in starting jet experiments in order 
to calculate the very important parameter of ambient 


air entrainment. For example, Cossali et al (2001) used 


Laser Doppler Anemometry to map the velocity field in 
a large region of a transient jet having Reynolds num¬ 
ber Re^ 12 * 10^ and Witze (1980) performed hot -film 


anemometer velocity measurements on the centreline 
of a starting jet at Re^ 7 * 10^ . Also, computational 
studies of gaseous starting jets have been performed 
for similar Reynolds numbers, e.g. Post et al (2000). 


In particular, Hu et al (2012) calculated the starting 


jet experiment by Witze (1980) in order to gain fur 


ther understanding of the behaviour of transient jets 
as related to internal combustion engines. However, the 
relevant Reynolds number in, for example, high pres¬ 
sure diesel spray injection depends on the prevailing 
conditions in the fuel nozzle and in the combustion 
chamber. Assumptions for characteristic quantities give 


a Reynolds number ^ 15 * 10^, e.g. jde Villiers et al 
(2004), or ^ 150 * 10^ when using state methodologies 
for the calculation of viscosity at supercritical condi¬ 
tions in modern engines, Oefelein et al ( |2012| ). 

The structure of the paper is as follows. Sectionj^de- 
scribes the flow field and the optical measurement tech¬ 
nique. Section [^describes the Wiener filter along with a 
calculation of the ensuing filtering errors and presents a 
comparison with another dissipation correction method 
Tanaka and Eaton ( 2QQ7| ). Section [^provides ; 


^ an assess¬ 


ment of the effect of filtering and demonstrates scalar 
dissipation rate measurements in a starting jet flow. Fi¬ 
nally, section gives a summary and conclusions. 


2 Experimental arrangement 

This section describes the flow configuration, the op¬ 
tical arrangement for laser induced fluorescence mea¬ 
surements and gives details for the measurement spatial 
resolution in relation to the flow length scales. The flow 
is that of an unsteady gaseous jet inside a low velocity 
CO- flow. A starting jet is created by suddenly releasing 
fluid from a nozzle. It has received relatively little at¬ 
tention compared with the steady state jet. However, it 
appears in a variety of important engineering applica¬ 
tions. For example, the injection inside the combustion 
chamber of a compression-ignition internal combustion 
engine is generally short, so that the fuel jet does not 
reach a steady state. In such engines, the mixing af¬ 
forded by the injection event is crucial in controlling 
emissions. 


2.1 Flow configuration 

The flow is established inside a square chamber of 30*30 
cm^ cross section. A general view of the flow chamber 
and of the coordinate system used is shown in figure 
The bottom part of the chamber serves as the co-flow 
conditioning area and the upper part, which is open to 
the atmosphere, is the injection and measurement area. 

The CO- flowing air, provided by a compressor, is 
heated using inline air process heaters, settles in a man¬ 
ifold and feeds the flow chamber through 8 inlets. The 
average cross-sectional temperature of the air inside the 
flow chamber is 65 °C±7 °C, in order to prevent acetone 
condensation, up to a downstream distance of 100 mm, 
equal to 25 nozzle diameters, from the jet nozzle exit. 
The temperature during data collection remained prac¬ 
tically constant and the above variation consisted of 
differences in temperature across the chamber, with 
the temperature near the chamber walls being slightly 
smaller. The jet injection and the optical measurements 
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Fig. 1 The experimental arrangement and the coordinate 
system used. 


take place after a flow conditioning section that ensures 
a low velocity (around 0.1 ™/s), low turbulence, homo¬ 
geneous velocity profile. Due to the low velocity of the 
co-flow, we estimate that air entrainment inside the in¬ 
jector nozzle, in-between injections, is minimal. 

The injector is a commercial, automotive gas injec¬ 
tor from KEIHIN (model KN3-2B). It consists of a fast 
acting solenoid valve, which is placed 20 mm upstream 
of the nozzle exit, and a straight tube of d=4mm in¬ 
ternal diameter, with a lip thickness of about 1 mm. 
It is connected to a settling chamber which, in turn, is 
connected to a pressurised air bottle through a pres¬ 
sure regulator that keeps the pressure inside the set¬ 
tling chamber constant. A pressure gauge and a K-type 
thermocouple are used to monitor the pressure and the 
temperature inside the settling chamber, upstream of 
the injector. The back pressure of the injection inside 
the settling chamber is 2.5 bar. 

The jet air flow passes, first, through inline process 
heaters and, later, through an acetone bubbler before 
being injected. The acetone bubbler is kept inside a wa¬ 
ter bath at 50 °C in order to maximise the acetone sig¬ 
nal and keep the level of seeded acetone constant dur¬ 
ing a run, whereas the temperature of the air stream 
downstream of the bubbler and just before the injector 
is kept at 65 °C. The co-flow air starts to flow well in 
advance of any measurements, so that thermal equilib¬ 
rium between the injector nozzle and the co-flow air is 
achieved. The concentration of acetone in the jet air is 
such that there is no discernible absorption of the laser 
light along the diameter of the jet, at distances very 
near the nozzle exit. 

In order to identify the temporal response of the in¬ 
jector’s solenoid valve opening, the acetone seeded fluid 
exiting the injector is imaged and the time delay be¬ 
tween the solenoid valve activation signal and the first 
camera signal is determined. Then, the ’’travel” time of 


the fluid in the straight pipe section between the valve 
and the nozzle exit is subtracted and the resulting time 
is considered as representing the solenoid valve opening 
time. The ensuing time is less than 1 ms. 

The injection velocity is determined from the cali¬ 
bration curve provided by the manufacturer of the in¬ 
jector and gives the output flow rate as a function of 
back pressure (pressure inside the settling chamber) 
and injection duration. In the absence of detailed ve¬ 
locity measurements, the injection velocity profile as 
a function of time (the injection schedule) is assumed 
to be a top hat profile, i.e. the fluid velocity is con¬ 
stant throughout the injection event and this velocity is 
reached instantaneously. The injection duration is 10 ms 
and the first appearance of jet fluid at the nozzle exit 
happens at about 5 ms after the initiation of the injec¬ 
tion pulse to the injector driver. 

The injection velocity is tx=20 ™/s and the Reynolds 
number, assuming a kinematic viscosity z^=1.8*10“^ 
is Re=4500 . From the velocity and the nozzle diameter 
a time scale is calculated as T = d/tx=0.2ms. 

During every injection, one image is acquired at a 
given time After the Start of Injection (ASI) and at 
every time ASI 300-500 images are acquired, in order 
to collect time-dependent, ensemble-averaged measure¬ 
ments. Injections are taking place every 500 ms, so as 
to avoid interaction between two consecutive injections 
and prevent acetone condensation inside the air lines. 
Images are acquired every 4th injection. 


2.2 Optical measurement 

The 4th harmonic of a pulsed, NdiYAG laser at 266 nm 
(10 ns pulse duration) is used to excite the acetone 
molecules and a 16-bit Intensified CCD camera (An- 
dor DH534-18F-03) is used to record the fluorescent 
intensity. A BC-3 (Schott) colour filter rejects light at 
the laser wavelength and allows the fluorescent light at 
the blue region of the spectrum to pass through. Three 
cylindrical fused silica lenses with focal lengths -20 mm, 
104 mm and 310 mm are used to generate a 45 mm high 
and 130 /im thick laser sheet. The camera is fitted with 
a 50 mm, f#1.2 (f^^ is the lens f-number) Nikkon lens 
and with a 13 mm extension ring for tighter focusing. 
Before entering the flow chamber, a beam-splitter re¬ 
flects 5% of the laser sheet energy into a dye cell con¬ 
taining a solution of Rhodamine 6C. A camera captures 
the emitted fluorescence signal, so that the laser energy 
and the laser sheet profile are monitored on a shot- 
by-shot basis. The laser energy at the probe volume is 
lOOmJ/pulse. A plan view of the optical setup is given 
in figure 
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16-bit 1024*1024 pixels^ 


Intensified CCD Cameras 


Fig. 2 Plan view of the optical setup. 


For conditions of constant temperature and pres¬ 
sure and at the optically thin limit, the acetone fluo¬ 
rescent intensity is proportional to, and varies linearly 
with, the laser energy and the acetone concentration. 
The images from the dye cell are used to correct the 
acquired acetone LIF images for the laser sheet pro¬ 
file and energy; typically, the inter-shot variation of the 
laser energy was of the order of 2% and, after these cor¬ 
rections, the spatial variation of the fluorescent signal 
along a radial profile at the exit of the nozzle was of 
the order of 2%; the camera non-linearities were min¬ 
imal and were not corrected for. The fluorescent sig¬ 
nal can further be processed to provide quantitative 
measurements of the mixture fraction. In the present 
two-fluid system the mixture fraction is the contribu¬ 
tion of the jet fluid stream to the local concentration. 
The corrected LIF signal intensity at any given spa¬ 
tial position in the flow is I{x,y) and the acetone LIF 
signal within the potential core at the nozzle exit is 
Iq. This intensity corresponds to pure jet fluid, so the 
mixture fraction is given at every point in the flow as 
C{x,y) = I{x^y)/Io^ where ( is the mixture fraction. In 
all images, the field of view includes the nozzle exit, 
so the acetone LIF signal corresponding to jet fluid 
is always known. The 2-dimensional approximation to 
the scalar dissipation rate is obtained from the mix¬ 
ture fraction field by using second order finite differ¬ 
ences to calculate the spatial derivatives, for example, 
d(/dx = (Ci+i,i ~ where Ax is the pixel 

size and Ci+i,i = C + ^^9 v )- fhe choice of the differ¬ 
entiation scheme depends on the resolution and signal- 
to-noise ratio of the measurements. While a higher or¬ 
der differentiation scheme has better spectral response, 
its outcome depends on the relative length scales be¬ 
tween the scalar structure and the spatial resolution. 
We, also, used a scheme that rotates the pixel array by 
45 degrees and calculates the derivative as an average of 
the original and rotated fields. The resulting difference 
in the mean scalar dissipation rate is less than 2% from 
the values obtained through 2nd order differencing. 


2.3 Spatial resolution and flow length scales 


The spatial resolution of the measurement depends on 
the laser sheet thickness and the in-plane spatial resolu¬ 
tion. The in-plane spatial resolution is not necessarily 
the same as the pixel spacing (which is 0.045 
for the 1024*1024 pixels^ CCD array) due to the blur¬ 
ring effect of the imaging system (camera lens, pixel 
size and intensifler). Using standard optics methods 
Siegmanj ( |1986 ), [Eckbreth (1996), the laser sheet thick¬ 
ness is estimated to be 130/im, at the 1/e^ position 
of the gaussian laser beam, growing to 150 ym at the 
edges of the imaged region - we also measured the laser 
beam thickness, using a scanning knife edge method, 
and found that it was ^150 ym. 

The blurring of the optical system is characterised 
by the Point Spread Function (PSF) (and its Fourier 
transform, the Modulation Tranfer Function - MTF) 


Goodman (1996), which quantifies the spreading of in¬ 
tensity values over the pixels when imaging a point 
source. The recorded intensity, of an imaged 

object, f{x, y), is given by fM{x, y) = y^h{x-x',y- 
y')f{x'^ y') dx'dy' ^ where h(x, y) is the PSF. In the case 
where the resolution of the optical system is limited by 
the camera lens, the PSF would be expected to have a 
form similar to gaussian. Where an image intensifier is 
employed, the form of the PSF might have longer tails 
than a gaussian. We measured the resolution of the op¬ 
tical system using the scanning knife edge technique. 
Figure shows the measured data points along with 
a spline fit and the derivative of the fit, the so-called 
line spread function. The Fourier transform of the line 
spread function is the MTF along a given direction, so 
it is closely related to the PSF. The resulting full width 
at half maximum of the LSF is 0.24 mm, however the 
LSF shows ’’tails” that are longer than an equivalent 
gaussian profile, which is sometimes used for fitting the 
LSF. 

The smallest length scale of variation of the velocity 
field in a turbulent flow is the Kolmogorov length scale 
defined as ^ where v is the kinematic vis¬ 

cosity of the fluid and 5 is the kinetic energy dissipation 
rate. 

With respect to a passive scalar field the smallest 
scale of motion depends on the relative magnitudes 
of momentum and mass diffusivity, as given by the 
Schmidt number. Sc = v/D {D the molecular dif¬ 
fusivity). The smallest length scale of the scalar field 
is then given by the Batchelor length scale as r]B = 


Batchelor 


(1959), 


Batchelor et al 


(1959) 


Corrsin 


(1964). When Sc^l, as in the present mixing of two 


incompressible gases, the Kolmogorov and Batchelor 
length scales are the same. 
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Fig. 3 Measurement of the optical system resolution. The 
symbols correspond to the measured intensity distribution 
at a pixel as a function of the knife edge position along the 
CCD array and the solid line is a spline fit to the data points. 
The derivative of the fit is the dashed line, known as the line 
spread function. 


-2 



Wavenumber [mm 


Fig. 4 Determination of the Batchelor length scale from the 
measured dissipation power spectra. The vertical arrows show 
the wavenumber corresponding to 2% of the peak value. 


A way to estimate the Batchelor length scale of the 
flow is to use the scalar energy spectrum Wang et al| 
(2007b). According to a model spectrum Pope (2000) 
the wavenumber, k, where kt] = 1 corresponds to a 
value of the dissipation spectrum equal to 2% of the 
peak value. The dissipation spectra are calculated in 
the following way: a small square region (around 50 
pixels on each side) at a position in the flow and at a 
given time ASI is chosen and the spectra of the mixture 
fraction axial gradient for all columns in the region, and 
for all images, are calculated and averaged together so 
as to obtain an ensemble-averaged spectrum of d(/dx 
in the axial direction; the same happens for all rows of 
the region, corresponding to the spectrum of dC^jdy in 
the radial direction. Before calculating the spectrum, 
a 2D Hann window is applied to the derivative data. 
This helps to minimise edge effects that arise from the 
abrupt cutoff of the signal as a result of the windowing 
of the data. Figure shows these spectra for a flow re¬ 
gion at the spatial location of the average position of the 
core of the vortex ring at time t/T=13.4 ASI and, also, 
shows the wavenumber corresponding to 2% of the peak 
value of the spectrum. For the axial-derivative spec¬ 
tra the wavenumber is ^ 3.3 mm“^ that corresponds 
to an estimated Batchelor scale of ^ 300 ym and for 
the radial-derivative spectra the corresponding num¬ 
bers are ^ 3.1 mm“^ and ^ 320/rm, respectively. For 
other times ASI and other flow positions the values 
range from ^ 250 ym to ^ 320 ym. So, taking as an es¬ 
timate of the Batchelor length scale the smaller of the 
two, r]B ^ 250 jam. So, these estimates confirm that the 
probe volume size offers appropriate spatial resolution 
for the measurements. 


Finally, the signal-to-noise ratio of the measurements 
can be deduced using the scalar energy power spectra 
calculated above. Exponential functions are fitted at 
the high wavenumber portion of the spectra which per¬ 
mit the calculation of the ratio of the integral of the 
fitted spectrum over the integral of the measured mi¬ 
nus the fitted spectra. These areas are shown in figure 
[^as the darker and lighter shaded areas, respectively. 
This ratio is the signal-to-noise ratio (SNR), which was 
always SNR>50 and typically SNR^^lOO, depending on 
the position in the flow and the time after the start of 
injection. Following the same procedure for the scalar 
gradient power spectra results in SNR^IO, i.e an or¬ 
der of magnitude lower than in the case of the mixture 
fraction. These numbers show that the mixture frac¬ 
tion measurements can be considered of good quality, 
however still the determination of the scalar dissipation 
rate is problematic. 


3 Digital filtering 

This section presents the theory behind the Wiener fil¬ 
ter and shows how to calculate the remaining errors, 
which can be considered as an estimate of the accuracy 
of the scalar dissipation rate measurements. Finally, the 
Wiener filter is compared with an analytically-derived 
method which corrects for noise and resolution effects 
and gives an estimate of the accuracy, in order to ver¬ 
ify the validity of the assumptions used in the Wiener 
filter. 
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Fig. 5 The spectrum of the scalar fluctuations for a small 
region in the flow centred at (x/D,y/D) = (4,1.15), showing 
the noise contribution (light fllled area) and the signal, with 
an estimated high-wavenumber variation (dark fllled area). 
The ratio of those two areas deflnes the signal-to-noise ratio 
of the measurements, at this particular location. 


3.1 Optimal Wiener filter 


In the following discussion, lowercase letters denote in¬ 
tensity images in physical space and uppercase letters 
the corresponding images in wavenumber space. The co¬ 
ordinates in physical space are r = (x^y) and the corre¬ 
sponding symbols in wavenumber space are k = {u,v). 

An image blurred by a Point Spread Function (PSF) 
and corrupted by additive noise can be written as 

-|-oo 

fM{r)= J h{r - r')f{r')dr'+ n{r) ( 2 ) 

— OO 

or Fm{i^) = H^N(k) , where is the mea¬ 

sured image, h{r) is the PSF of the measuring system, 
/(r) is the noise-free, unblurred image, and n{r) is the 
zero-mean noise. Using the above measurement model, 
the Wiener filter tries to reconstruct the ideal image, 
/(r), by applying a filter g{r) to the measured image 
Then, the filtered image can be written as 

-|-CX3 

frir) = j 9{r- r')fM{r') dr' (3) 

— C» 


or Ff{iF) = G{k)Fm{i^)- The filtering error is defined 
as e(r) = f{r) — /f(^) and the mean square error is 


= E 




( 4 ) 


where E denotes the expectation. The Wiener filter is 
the optimum solution that minimises this mean square 
error. A complete exposition of this minimisation prob¬ 
lem is given in Petrou and Bosdogianni (1999), where 


it is proven that the Fourier transform of the function 
g{r) is G{hz) = W{n)/H{n)^ where 


W{n) = 






( 5 ) 


is the Fourier transform of the Wiener filter. The power 
spectral densities of the noise and of the uncorrupted, 
noise-free image are 5'^(/i:) = |A/'(/i:)p and 5'(/^) = \F{k)\‘^, 
respectively and the filtered image is 


Ff{^) 


W{k)Fm{^) 

H{k) 


(6) 


A characteristic of the Wiener filter is that it re¬ 
quires knowledge of the power spectral densities of both 
the unblurred, noise-free image and the noise. A possi¬ 


bility, followed by Krawczynski et al (2006a) in the con¬ 


text of fluid mechanics and, also, suggested by [Petrou 


and Bosdogianni (1999) for a more general situation, is 


the following: a series of measured images can be used 
to plot an estimate of the spectral density of the sig¬ 
nal plus noise. With increasing wavenumbers the power 
spectral density of the signal decreases, up to the point 
where the noise starts to have a significant contribution. 
Then, the spectrum of the true signal, for the highest 
wavenumbers, can be modelled by a decreasing expo¬ 
nential function. Such a high-wavenumber form of the 
three dimensional spectrum has been proposed at var¬ 


ious model spectra, Corrsin (1964); Pao (1965); Batch¬ 


elor et al (1959); Pope (2000). Similarly, a suitable fit 


{e.g. using cubic splines) can be found for the noise den¬ 
sity and extrapolated to the wavenumbers where the 
true signal dominates the power spectrum. 

Figure shows an ensemble-averaged 2D spectrum 
of the mixture fraction fluctuations at time t/T=9.65. 

At the low wavenumbers the contours of the 2D spectra 
are not circular which demonstrates the in-homogeneity 
of the flow field at the large scales of motion. Never¬ 
theless, at higher wavenumbers the contours become 
circular and the exponential fitting only takes place 
after this radial symmetry has been established. The 
models used for the construction of the Wiener Alter 
at t/T=9.65 are shown in figure which plots a hor¬ 
izontal cut through the centre of the frequency plane 
and along the radial direction of the averaged 2D spec¬ 
trum of figure alongside the model spectra of the 
high wavenumber region of the signal and of the noise. 

The portion of the spectrum that is modelled in the 
determination of the Wiener Alter will be discussed in 
section 14.11 

The Wiener Alter calculated from Eqj^at time t/T=9.65 
is shown in figure and it shows a ” bump”, as a char¬ 
acteristic of the Wiener Alter with the inclusion of the 
PSF. 
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Fig. 6 An averaged 2D spectrum of the mixture fraction 
fluctuations at time t/T=9.65. 



Fig. 8 The transfer function of the Wiener Alter at time 
t/T=9.65. 


on k) 


8 



Wavenumber [mm 


Fig. 7 A horizontal cut through the averaged 2D spectrum 
of the mixture fraction fluctuations at time t/T=9.65. The 
model signal spectrum for the high wavenumbers, the model 
spectrum for the noise and the modulation transfer function 
(which is the Fourier transform of the PSF) are also shown. 


+ 00 

e" = y A [i7|(vr - 1)2 + N^w^] dK (8) 

— OO 

where Fs = HF. If the optimum filter from Eq. is 
substituted in the above equation, the mean square er¬ 
ror becomes 


-\-oo 


* dK 


F| + 


( 9 ) 


The mean square error in the spatial derivative is cal¬ 
culated by considering the following mean square error 
form 


+ OOr 

r 


dfp 

2 

1 


dfp 

2~ 

J 

-OO 

dx 

dx 


dy 

dy 



( 10 ) 


3.2 Filtering errors 


where /dx — d/p/Oxf dr and 


The derivation of the Wiener filter provides the pos¬ 
sibility to calculate the ensuing filtering errors for the 
variance of the mixture fraction and the scalar dissipa¬ 
tion fields. By using Parseval’s theorem, a different way 
to express the mean square error, Eq. is 

-\-oo +00 

e2= y |/(r)-+(r)|2dr= J |F(ac)(k)| 2 dAc (7) 

— OO —OO 

Using Eqs. and and taking into account that the 
signal and the noise are uncorrelated, the mean square 
error, Eq.[^ can be written as (omitting the dependence 


+ 00 +00 

4, = / = / .^ |FW - 

— OO —OO 

and similarly for e‘^ y. So, 

+ 00 

= J \F{k,) — FF{K,)f dK, 

— OO 

where k? = v‘^. The expression \F{k) — FF{K)f 

appears in the mean square error of Eq.j^and is given in 
Eq.[^ when the optimal Wiener filter is used. Then, the 
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associated mean square error for the spatial derivative 
of the mixture fraction field is 


+ 00 


el = 


1 






dK 


( 11 ) 


The mean square filtering errors were evaluated from 
Eqs. and pT] and are shown in table where ( is the 
instantaneous mixture fraction, for all the times after 
the start of injection. The mean square errors of the 
mixture fraction field are given as the ratio of the Eq. 

and the integral over all wavenumbers of the model 
spectrum, d/i:, i.e. the ’’true” variance. The 

corresponding mean square errors of the scalar spatial 


gradient are given as the ratio of Eq. and the inte¬ 
gral over all wavenumbers of the model dissipation spec¬ 
trum, dK. Erom Eq. the mean square er¬ 

ror can be thought of as an excess variance contributing 
to the measurement of the true variance. Similarly, the 
mean square error of the derivative, Eq.[^ is the corre¬ 
sponding excess variance to the variance of the spatial 
derivative, i.e. the scalar dissipation. In summary, ta¬ 
ble shows that average error of the mixture fraction 
variance over all times ASI is less than 2% and that the 
scalar dissipation rate is accurate to around 20%. This 
level of accuracy can be compared to values reported in 
the literature. Eor example, in a confined jet experiment 
the budget of the scalar variance gave an uncertainty 


within 20% Markides and Mastorakos (2006) and in a 


partially stirred reactor Krawczynski et al (2006b) the 


accuracy was again within 20%, judged by the scatter of 
the values by using different methods to calculate the 
dissipation rate. The possibility to estimate explicitly 
the accuracy, Eq. IE and the sensitivity of the filter¬ 
ing process, Eq. below, without relying in ad-hoc 
assumptions render the present technique suitable for 
post processing scalar measurements with the aim of 
calculating the scalar dissipation rate. 

Einally, we comment on the effect of the spectral 
fittings that are used to define models for the noise 
spectrum and the high wavenumber signal spectrum 
and show that suboptimal models affect little the de¬ 
rived quantities. Suppose that the Wiener filter was es¬ 
timated with model spectra that were different from the 
real spectra, thus producing a filter function IT* instead 
of the optimum one, W. Then, the mean square error 
(ej) for this choice of Wiener filter is given by Eq. 
where IT is replaced by IT*. Expanding the expression 
inside the braces and then completing the square in the 
ensuing expression results in 


el = 


+ 00 

r 


+ (IT*-IT)V 


1 




F2 


dK 


The last term inside the braces in the above equation 
is the mean square error for the optimum Wiener filter, 
Eq. So the mean square error for the ’’sub-optimal” 
Wiener filter is given by 

+ 00 

el = e^+ f (w, _ wf 

— OO 

This expression is quadratic in the difference between 
the two Wiener filters, implying that even a non-perfect 
evaluation of the model spectra is still capable of pro¬ 
ducing small filtering errors, comparable with the min¬ 
imum one. Eor example, IT* = (1 -h (a)IT, where a is 
assumed constant for all wavenumbers. Erom Eq. it 
follows that W^/H^ = Fl/Fl^ = F|/ (F| + N^) and 
the above equation becomes 

+ 00 

el = e^ + a^ j F^df ( 12 ) 

— CX) 

SO the extra mean square error for, say, a = 0.1 (or 
10%) is 1% of the measured variance for both the mix¬ 
ture fraction and the scalar spatial derivative. There¬ 
fore, the current analysis shows that the Wiener filter 
is rather insensitive in the choice of models, a fact that 
was verified for some cases by choosing alternative mod¬ 
els for either the measurement or the noise spectra. 

In summary, the Wiener filter tries to recreate the 
noise-contaminated turbulent signal by forcing an agree¬ 
ment (in the least mean square sense) between the mea¬ 
sured and the ’’ideal” spectra. In doing so, the scale- 
dependent effect of noise is taken into account in order 
to reconstruct the noise-free signal. It is still important, 
however, to analyse the effectiveness of the process and 
this is done in section [4j 


3.3 A dissipation correction method based on 
Richardson’s extrapolation 


The method of Tanaka and Eaton (2007) provides a cor¬ 
rection to energy dissipation rate measurements (easily 
extendible to scalar dissipation rate as well) for the ef¬ 
fects of noise and inadequate spatial resolution. Since it 
can be shown that this method is second order accurate 
in the spatial resolution (thus giving an estimate of the 
accuracy), we briefly present it here and use it in sub¬ 
sequent sections to compare it with the Wiener filter, 
in trying to render greater credibility to the modelling 
assumptions of the Wiener filter adopted earlier. This 
method, however, can only correct for the mean dissipa¬ 
tion rate, so it cannot be used when the instantaneous 
fields are required. 
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Table 1 Percentage (%) errors of the variances of the mixture fraction and of the scalar spatial derivative, after using the 
Wiener filter. The errors are normalised with the model scalar variance and the model scalar dissipation, respectively. 


Time ASI, t/T 

0.9 

2.15 

3.4 

4.65 

5.9 

7.15 

9.65 

10.9 

13.4 

15.9 

(C^> 

3.5 

1.9 

1.6 

2.4 

2.3 

0.8 

2.0 

1.0 

1.5 

1.3 

((VC)"> 

19.2 

16.1 

29.2 

18.2 

23.3 

10.7 

19.9 

13.9 

22.6 

23.6 


The idea of Richardson extrapolation is that by re¬ 
peating a numerical calculation with different values of 
a parameter (in this case the spatial resolution), the 
’true’ value of the quantity that needs to be calculated 
can be obtained by extrapolation to zero value of the 
parameter. In Tanaka and Eaton (2007) this approach 
is used to calculate the turbulent kinetic energy dis¬ 
sipation rate from Particle Image Velocimetry (PIV) 
images, where the spatial resolution is adequate and 
experimental noise dominates. The energy dissipation 
rate is evaluated at two spatial resolutions. Ax and 
2Z\x, and a combination of these two measurements 
gives a measurement of the ensemble averaged dissipa¬ 
tion rate which is second order accurate on Ax. It is 
straightforward to adapt this correction method to the 
scalar dissipation rate and the formula to calculate the 
mean scalar dissipation rate is 


ix) = iXA) 

where X 2 /i = {^CI^x)l^^ + {A(jAy)l^y and the terms 
are calculated as, for example, 


f \ _ Ci+2,j Ci-2,j 

\^^J2Aco 4Zla; 

where i^j denote the pixels and Ax, Ay the pixel spac- 
ings in the two directions. 


4 Results 

In this section the effects of filtering on the measured 
mixture fraction distributions are presented and results 
from the starting jet flow are reported in order to demon¬ 
strate the measurement approach. 


4.1 Filtering model noise 

The portion of the spectrum that is modelled through 
the exponential decay in the construction of the Wiener 
Alter should be as small as possible. Furthermore, the 
properties of the residual fields (measured minus Al¬ 
tered) have to be tested to check the effectiveness of 
the Altering operation. 

The level of modelling of the measured spectrum 
(Agure[^, when determining the Wiener Alter, can be 


given by the portion of the measured variance that is 
modelled by the large-wavenumber exponential decay. 
Using the same procedure as earlier (c./. Agure|^ this 
portion is given by the ratio of the area between the 
measured and Atted spectra over the area under the 
measured spectrum. The ’cut-oA’ wavenumber, where 
the modelled spectrum starts to depart from the mea¬ 
sured spectrum is of the same order as the smallest es¬ 
timated Batchelor length scale, which is true for other 
times ASI as well, showing that the modelling does 
not aflect the smallest scales of the scalar held. Table 
shows the relevant percentages of modelled variance 
and the ’cut-oA’ wavenumbers, for all times ASI. The 
mixture fraction variance is almost entirely determined 
by the measurements, with the Atted part of the spec¬ 
trum being less than 1% in all cases. This should be 
expected given the relatively high SNR calculated ear¬ 
lier. The large ’cut-oA’ length scales for the Arst couple 
of times ASI are a consequence of the development of 
the scalar held at these earlier times, where the fluid 
has just left the nozzle and the scalar held has not had 
time to develop turbulent fluctuations. The spectrum 
is determined only by the very large scales, with no 
intermediate scales present, so the noise contribution 
extends at relatively low wavenumbers. 

In assessing the statistical properties of the correc¬ 
tions made to the raw images, Agure[^ presents a raw in¬ 
stantaneous mixture fraction image at time t/T=9.65, 
the same instantaneous image after the application of 
the Wiener Alter and the diAerence between these two 
images (raw minus Altered). 

The two instantaneous images are practically the 
same to each other (as shown later the correction is 
less than 2% of the measured mixture fraction values 
and, typically, less than 1%) and the general charac¬ 
teristics of the evolving jet are not changed by the ap¬ 
plication of the Alter. In the Altered image the contour 
lines are smoother without altering their location and 
shape. Qualitatively, this shows that the Wiener Alter 
is adequate in removing noise contributions, with more 
quantitative results given later in this and the next sub¬ 
sections. 

Initially, the noise characteristics of the imaging sys¬ 
tem are assessed from the recorded images. The view¬ 
ing area in all images contains the jet plus regions of 
the laser sheet where no acetone vapour is present; in 
these regions the acquired signal should be constant 
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Table 2 The cut-off wavenumber (mm ^), where the model spectrum starts to deviate from the measured spectrum and 
percentage (%) of the measured spectrum of the mixture fraction that is modelled in the Wiener filter. 


Time ASI, t/T 

0.9 

2.15 

3.4 

4.65 

5.9 

7.15 

9.65 

10.9 

13.4 

15.9 

^cut-off 

1.2 

1.9 

3.5 

2.5 

3.0 

2.7 

2.9 

2.8 

1.7 

2.6 

-F(k.) modelled 

0.26 

0.08 

0.01 

0.03 

0.02 

0.01 

0.01 

0.01 

0.03 

0.02 




y/D 


Fig. 9 An instantaneous raw mixture fraction image at time t/T=9.65 (left), the corresponding instantaneous image after 
application of the Wiener filter (centre) and the difference between the instantaneous images, filtered minus raw (right). The 
scale corresponds to mixture fraction. 


with only the camera noise superimposed. For exam¬ 
ple, subtracting the ensemble averaged image from a 
raw instantaneous image and collecting the pixel val¬ 
ues inside the rectangle shown in the leftmost plot of 
figure will give the noise distribution of the camera. 
Figure plO] shows the distribution of the noise samples 
from this region compared with a standard normal dis¬ 
tribution (the quantile plot is constructed by plotting 
ordered data values of two distributions against each 
other; in case the two data sets follow the same distri¬ 
bution the plot is a straight line). For most part of the 
graph the noise sample appears to be close to the nor¬ 
mal distribution, at least up to 2 standard deviations 
and noise samples from different instantaneous images 
and for different times ASI show a similar behaviour. 
So, we can assume that the noise characteristics of the 
camera follow a gaussian distribution. 

The image obtained as a result of the filtering op¬ 
eration is the best estimate of the unblurred, noise-free 
image. So, the noise field (the difference between the fil¬ 
tered and the raw images - c.f. the image model, Eq§ 
should be both as small as possible and as close to white 
noise as possible. It should be arguably accepted that 
the level of the correction should be as small as pos¬ 
sible and in line with the signal-to-noise ratio. Also, 
the filtering should preserve the features of the mea¬ 
sured flow field and not remove any large or smaller 
scale characteristics. So, it should be a featureless noise 
field, i.e. white noise. The difference image in figure 
1^ shows some remaining ’structure’ after filtering. The 


-3 



Fig. 10 Quantile-quantile plot for the camera noise. The 
noise sample is produced from an image of the deviations 
from the ensemble average for a region outside the jet flow. 

level of this residual structure is small, less than 3% 
of the measured value, however it points, as explained 
later, to one of the limitations of the filtering technique. 

A noise distribution is shown in figure [TT| The sam¬ 
ple is obtained by collecting all the pixels inside the 
small rectangle in the rightmost image of figure the 
difference between filtered and measured images. The 
noise distribution compares well with the normal distri¬ 
bution with only minimal deviations from the gaussian 
distribution at fluctuations larger than about 2 stan¬ 
dard deviations. The distributions at other positions in 
the flow (from regions within the jet contour) and at 
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other times ASI are very similar. This particular noise 
distribution has a mean value of zero and variance 10“^; 
in general, the level of the correction is less than 1% of 
the measured value of the mixture fraction. This level 
can be assessed by considering, for example, that when 
the noise is gaussian and additive, at every spatial lo¬ 
cation in the image its distribution is /SNR)^ 

where N {a ^b) is the normal distribution with mean a 
and variance b and a is the standard deviation of the 
measurement at the same spatial location. So, for a typ¬ 
ical value of the standard deviation of the mixture frac¬ 
tion fluctuations a ~ 0.1 and SNR ^ 100 the variance 
of the noise should be ~ 10“^, as measured. Conse¬ 
quently, the correction imposed by the Wiener filter is 
generally small and consistent with the signal-to-noise 
ratio of the measurements and the distribution of the 
residual noise samples follows well a normal distribu¬ 
tion. This means that the Wiener filter does not alter 
the underlying mixture fraction field and, more impor¬ 
tantly, does not introduce any bias in the values of the 
corrected images. This figure also verifies the assump¬ 
tion that the camera noise follows a gaussian distribu¬ 
tion. 



Fig. 11 Quantile-quantile plots of the optimal Wiener filter 
model noise compared with a standard normal distribution. 


The 2D spectrum of the noise values in the same 
small rectangle of the rightmost image of figure is 
shown in figure This spectrum has been calculated 
from a single image, which explains the random fluctu¬ 
ations evident in this plot and the ensemble averaged 
spectrum shows the same characteristics as the single¬ 
shot spectrum. Small regions at other positions inside 
the jet body and for other times ASI give similar results. 
This spectrum is symmetric around the origin and at 
the larger wavenumbers the magnitude of the spectrum 
is almost the same everywhere, irrespective of orienta¬ 
tion or scale. This implies that there is no bias in the 


Altering of the very small scales. At larger scales, larger 
than about 0.4-0.5 mm, there is a small increase in the 
level of the spectrum. This scale corresponds roughly 
to the cut-off wavenumber where the model spectrum 
starts to deviate from the measured spectrum of the 
mixture fraction fluctuations, table At wavenumbers 
smaller than the cut-off wavenumber, the noise level is 
considerably smaller than the signal. However, the Al¬ 
tering operation still attenuates the measured signal. In 
this way, some signal structure is smoothed out, as seen 
in figure In contrast, at larger wavenumbers where 
noise starts to become comparable to the measured sig¬ 
nal the Wiener Alter acts to reconstruct the ideal signal 
by only removing the noise contribution. 



Radial wavenumber [mm ] 


Fig. 12 The 2D spectrum of the residual noise inside the 
rectangle in the rightmost image of figure The logarithm 
base-10 of the spectrum is shown in this figure. 


Finally, we note that the above discussion referred 
to flow regions away from the jet boundaries. Outside 
and inside the jet body the residual noise distribution 
is featureless and the ensuing corrections are small (as 
shown by the noise distribution in figure and the 
2D spectrum in figure . Despite that, the rightmost 
plot of figure shows that there is some dampening 
of characteristic features of the mixture fraction held 
at the jet boundaries. Since the jet boundaries always 
have very large gradients, these gradients will always 
be attenuated after applying the Alter. So, given the 
persistence of these gradient from injection to injection 
(in contrast to the random emergence of large gradients 
inside the jet body), the scalar dissipation rate will have 
a negative bias at the jet boundaries. 
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4.2 Scalar dissipation rate measurements 


This section presents results of the scalar dissipation 
rate measurements, compares the raw with the filtered 
images and proves the effectiveness of the Wiener filter 
in obtaining accurate scalar dissipation rate measure¬ 
ments. 

The scalar dissipation rate will be presented without 
the diffusivity, so its dimensions are in mm“^. We can 
assume that the diffusivity is constant throughout the 
flow with a value between 20mm^/s and 25mm^/s. In 
relation to figurej^ figure [T^ shows radial profiles, at the 
same axial distance, of both the raw and filtered image 
and, also, of the radial derivative of both signals. The 
Wiener filter proves effective in restoring the signal and 
removing small fluctuations that could be attributed to 
noise. 


In order to assess the ability of the Wiener filter to 
produce accurate results for the mean scalar dissipa¬ 
tion rate, figure shows the mean scalar dissipation 
rate at time t/T=9.65 for the Wiener filter and the 
correction method from Tanaka and Eaton (2007), as 
presented in section |3.3[ The two distributions are very 
similar to each other with differences being found in 
some structures within the jet body, which have very 
small dissipation values; it is, also, important to note 
that the values of the mean scalar dissipation rate, as 
calculated from the two methods, are almost the same. 
The Wiener filter practically, and as applied here, cor¬ 
rects for the noise in the images, whereas the correc¬ 
tion method from Tanaka and Eaton (2007) corrects 
for both noise and resolution effects. By taking into ac¬ 
count figure where the two methods give almost the 
same results, we can conclude that the spatial resolution 
of the measurements is enough to resolve the smallest 
scales of variation of the scalar field with a high signal- 
to-noise ratio and that we can, consequently, measure 
the scalar dissipation rate with good accuracy. 

Instantaneous scalar dissipation images of both the 
raw and the filtered measurements are shown in the 
top row of figure These distributions show a thin 
layer of very high scalar dissipation values surrounding 
the jet. Inside the jet body dissipation values are found 
in small regions within an otherwise calm fluid and the 
level of the dissipation values are somehow smaller than 
the values found at the jet boundary. The distribution 
of the scalar dissipation rate within the jet boundaries 
changes from injection to injection with the result be¬ 
ing a distribution of mean scalar dissipation rate having 
two peaks at jet shear layer with very low values in be¬ 
tween, as shown in figure Part of the thickness of 
the ensemble average scalar dissipation profile along the 
radial direction, especially for further axial distances. 


could result from large scale variations between injec¬ 
tions. We will quantify the effect, in subsequent analy¬ 
sis, by considering the repeatability of consecutive in¬ 
jection events through the distributions of the mixture 
fraction at locations near the nozzle exit, where interac¬ 
tion with the surrounding fluid is minimal; in general, 
the distributions at these locations are single-mode and 
compact. 


Qualitatively in the raw instantaneous image any 
structure is always hidden behind a noise Veil’. In con¬ 
trast, in the filtered image, this noise has been largely 
removed and the expected smooth regions of very low 
dissipation values appear embedded within larger and 
spatially concentrated dissipation values. Eor compar¬ 
ison, figure shows the corresponding filtered image 
assuming a delta function as the PSE. Including the 
measured PSE in the Wiener filter restores the mix¬ 
ture fraction gradients to larger magnitude, as opposed 
to assuming perfect imaging. Noise evident in the lat¬ 
ter image has, additionally, been removed as a result of 
the different form of the Wiener filter. The difference 
between the filtered and the raw scalar dissipation dis¬ 
tribution is, also, given in figure The largest correc¬ 
tion applied to the measured image in order to arrive at 
the filtered result is of the order of 40%. We can assess 
this level of correction by considering the apparent dis¬ 
sipation, e.g. Mi and Antonia (1994), which quantifies 
the effect of noise in the true scalar dissipation value. 

This quantity is defined as {xa) = 3 ^ = (Xm) - (xt) 
(considering, e.g., forward differencing for the deriva¬ 
tive), where subscripts a, m, t correspond to the appar¬ 
ent, measured and true dissipation values, respectively, 
is the noise variance and Ax is the spatial discretisa¬ 
tion. A rough estimate of this quantity for the present 
experimental setup results in ~ 30% of the measured 
value, which is comparable to the actual correction im¬ 
posed on the measurements. The differences between 
the measured and corrected images are mainly concen¬ 
trated at the edges of the jet and they are typically 
positive, meaning that the raw image gradients have 
been restored by the filtering operation; within the jet 
body the corrections have a negative bias. This be¬ 
haviour of the filtering operation results since the PSE 
acts to restore the large gradients, concentrated at the 
jet boundaries, whereas the attenuated gradients within 
the jet body are expected for a filtering technique like 
the Wiener filter, which treats the measurements in a 
’global’ manner, without being able to adapt to local 
conditions of the flow field. 
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Fig. 13 Radial profiles of an instantaneous raw mixture fraction image at time t/T=10.9 and x/D=4 and of the corresponding 
filtered image (left). Corresponding radial derivatives of the same profiles (right). 
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Fig. 14 T he mea n scalar dissipation rate at time t/T=9.65 for the Wiener filter (left) and the correction method from 
[and Eaton] ( |2007| (right). The images show the logarithm base-10 of the scalar dissipation without the diffusivity. 
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4.3 Discussion 

For any given experiment, the outcome of applying of 
the Wiener filter for the post processing of scalar mea¬ 
surements depends in an interconnected way on both 
the signal to noise ratio and the measurement spatial 
resolution. This is more clearly seen by rewriting the 
quantity W{n)/H{n) appearing in Eq. as 

^ ~H{^) l + Sn{K)/Ss{K) 

where Sg = H‘^S is the blurred, but noise-free, spec¬ 
trum (essentially the model spectrum in figure and 
the ratio Sn/Sg is the inverse of the signal to noise ra¬ 
tio as a function of the wavenumber. At the highest 
wavenumbers we expect the signal to noise ratio to fall 
off exponentially fast, and faster than the PSF, so the 
very large, and noisier, wavenumbers of the signal will 
be attenuated. At lower wavenumbers, the form of the 


ratio W{k)/H{k) depends on the relative fall-off of the 
’characteristic’ bandwidths of the PSF and the mea¬ 
surement. For example, the characteristic wavenumbers 
can be taken for PSF as being where the PSF starts to 
decrease and for the signal the wavenumber where the 
noise spectrum becomes important and thus where the 
measured spectrum diverts from the real one. If the real 
spectrum falls off earlier than the PSF, then the Wiener 
filter will roll off smoothly and will decrease the effect 
of the noise only. In the case when the bandwidth of the 
flow scales is comparable to the PSF resolution band¬ 
width, the Wiener filter will take a form as to initially 
accentuate the blurred image, before subsequently at¬ 
tenuating the noise components; in this situation some 
imaging artifacts might be created. In the case when the 
measurement resolution is very coarse, a lot of informa¬ 
tion will have been lost and the Wiener filter might not 
be the best solution to reconstruct the measured signal. 
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Fig. 15 An instantaneous raw scalar dissipation image at time t/T=9.65 (top left), the corresponding image after application 
of the Wiener filter (top right) and the difference between the filtered and the raw images (bottom left). Also shown is the 
instantaneous image assuming the PSF is a delta function (bottom right). The images show the logarithm base-10 of the scalar 
dissipation without the diffusivity; the difference image shows the difference of the dissipation values without the diffusivity. 


The qualitative discussion above is now made quan¬ 
titative by considering the filtering process for flows ex¬ 
hibiting higher levels of turbulence, by varying the rel¬ 
ative magnitudes of the flow length scales as compared 
to the experimental spatial resolution. We use the mea¬ 
sured PSF so as to retain the effect of the intensifler. 


We adapt a model for the 3D spectrum of veloc¬ 
ity Pope (2000) is adapted for use as a model for the 
3D spectrum of a scalar, as in Wang et al| ( |2007a ), 
and we obtain t he model ID spectrum as S(- (/^i) = 
(/^) Tennekes and Lumley ( 1972|, where 
E(^ (k,) is the 3D scalar spectrum. This rnodel ID spec¬ 
trum, the measured PSF and white noise are used in Eq. 
1^ to obtain the ’measured’ (blurred and noisy) spec¬ 
trum, having signal-to-noise ratio ^100 (to calculate 
the SNR the spectra are integrated up to the wavenum¬ 
ber corresponding to the Batchelor length scale). In or¬ 
der to model low and high turbulent conditions, the 
Taylor-scale Reynolds number takes the values 100 and 
1000. These quantities determine the parameters of the 
Wiener Alter using Eq. and allow the calculation of 
the Altering errors for the scalar and its dissipation us¬ 
ing Eqs. and IE respectively. As an example, figure 
16 , which is equivalent to figureshows the ID mod¬ 


eled, ’measured’ and extrapolated spectra along with 


the measured PSE for Taylor-scale Reynolds number 
1000. The resulting Wiener Alter is similar to the one 
presented in Agure 


6 



Fig. 16 The model ID spectrum, its blurred version, its 
blurred and noisy version and the measured PSF. The Taylor- 
scale Reynolds number is 1000 and — 1-25. 


For normalisation we use the value where the MTF 
is at 2% of the peak value, tjmtf = 0.24 mm, and in 
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figure r]B /vmtf = 1-25, which corresponds to the 
measured values of the MTF and flow length scales. 
In order, however, to check the effect of the separation 
between the fast roll-off of the measured MTF at high 
wavenumbers and the Batchelor length scale in addition 
to the experimental noise, we also use r]B /vmtf = 0.4 
and a range of signal-to-noise ratios. The filtering errors 
for all cases are given in table and figure 

In the cases in table [3] the characteristic scale PSF is 
either larger or smaller than the Batchelor length scale. 
As r]B becomes larger than the error in the scalar 

dissipation rate decreases, for both small and large Tay¬ 
lor Reynolds number. So, the first conclusion is that, 
for flow fields described by these spectra, as long as the 
Batchelor length scale is well resolved by the PSF, the 
Wiener Alter is capable of adequately Altering the mea¬ 
surements to produce a value for the scalar dissipation 
rate within acceptable limits. 

In an imaging experiment, increasing the Taylor- 
scale Reynolds number results in a reduction in the 
size of the smallest flow length scales and, if one keeps 
the experimental setup unchanged, both resolution and 
noise are negatively affected; so, the measurement er¬ 
rors become larger. The ideal case is a simultaneous in¬ 
crease in spatial resolution (increasing the experiment 
size is a potential avenue to achieve that, wherever pos¬ 
sible) and a proportional increase in signal power, per¬ 
haps by reducing the laser sheet size (even though this 
would eventually restrict the lowest resolvable wavenum¬ 
ber) . Alternatively, one can attempt to tweak either the 
spatial resolution or the signal level and a common ex¬ 
perimental situation, as the Taylor-scale Reynolds num¬ 
ber increases, is to keep the PSF and the SNR constant. 
This situation is presented in table in rows 1 and 4, 
where the spatial resolution reduces from tjb/Vmtf = 
1.25 to Vb/vmtf = 0.4 as the Taylor-scale Reynolds 
number increases from 100 to 1000. The resulting error 
in the estimate of the variance of the scalar decreases 
from 10“^% to 10“^%, whereas the error in the estimate 
of the scalar dissipation rate increases from 7.2% to 
30.6%. The decrease in the scalar variance error results 
from capturing a larger part of the low-wave number 
region of the spectrum as the Batchelor length scale 
decreases, while the imaging held of view remains the 
same. In contrast, the increase of the error in the scalar 
dissipation is expected, given the lower spatial resolu¬ 
tion. 

Although the common experimental situation is as 
described above, figure shows that it is more effec¬ 
tive to expand efforts in increasing the spatial resolution 
rather than to improve the SNR level. Figure pT| shows 
the effect of the signal-to-noise ratio at Rex = 10^ 
and 10^ when Vb/vmtf = 0.4 or 1.25. As expected. 


the error reduces with increasing signal-to-noise ratio; 
increasing the signal-to-noise ratio by order of mag¬ 
nitude decreases the error by a third in the case of 
Vb/vmtf = 0.4, whereas the error decreases by half 
when Vb/vmtf = 1.25. At the same time no increase 
in the signal-to-noise ratio seems to be able to achieve 
the error afforded when the flow length scales are better 
resolved; e.g. in figure pT| the line with Vb/vmtf = 0.4 
and Rex = 10^ (Ailed squares) is always higher than 
the line with Vb/vmtf = 1-25 at the same Rex = 10^ 
(Ailed circles). The implication is that attempting to 
achieve better measurement resolution at a given signal- 
to-noise ratio is the more profitable path to decreasing 
the error in the estimate of the dissipation. Further ev¬ 
idence of this point of view is provided by the com¬ 
parison of the error variation for Rex = 10^ between 
Vb/vmtf = 0.4 and Vb/vmtf = 1-25, in figure pT| Hav¬ 
ing a better resolution can counteract even an increase 
of an order of magnitude of the Taylor scale Reynolds 
number, irrespective of the signal-to-noise ratio. 



Fig. 17 The estimate of the error in the dissipation for differ¬ 
ent signal-to-noise ratios, for Taylor-scale Reynolds number 
Rex = 1000 and for vb/vmtf = 0.4 and vb/vmtf = 1.25. 


The preceding results corroborate the earlier qual¬ 
itative discussion that the relative magnitude of the 
measurement resolution cut-off with respect to the flow 
scales {e.g. the ratio Vb/vmtf) is the determining fac¬ 
tor in the successful application of the Wiener Alter. So, 
for the scales in this experiment, the spatial resolution 
is good enough to resolve the scalar dissipation rate us¬ 
ing the Wiener Alter, but, for a flow with smaller length 
scales an imaging system with a narrower PSF should 
be used. However, the Reynolds number by itself should 
not be the determining factor in deciding the use of the 
Wiener Alter in a particular situation. 
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Table 3 The filtering errors in the estimate of the scalar and its dissipation for a blurred and noisy model ID spectrum, with 
parameters the Taylor-scale Reynolds numbers, Re\^ and the Batchelor length scale, for signal-to-noise ratio 100. 


Rex 


SNR 

e2[%] 

4[%] 

100 

1.25 

100 

0.01 

7.2 

100 

0.4 

100 

0.05 

25.9 

1000 

1.25 

100 

5E-4 

12.7 

1000 

0.4 

100 

lE-3 

30.6 


5 Conclusions 


An implementation of the Wiener filter for the measure¬ 
ment of mixture fraction in a starting jet flow which 
is in-homogeneous at the largest scales of motion is 
presented. An explicit calculation of the accuracy of 
the mixture fraction variance and of the scalar dissi¬ 
pation rate showed that, in the present measurements, 
the Wiener filter produced results that are, on average, 
within 2% of the mixture fraction variance and within 
20% of the scalar dissipation rate. This accuracy es¬ 
timate was, also, checked using an independent correc¬ 
tion method and the results reinforced the applicability 
of the Wiener filter approach in the present measure¬ 
ments. We choose the Wiener filter because the error 
calculation gives an explicit estimate of the accuracy 
of the measurement of the scalar dissipation rate and 
because the Wiener filter can also correct the instanta¬ 
neous images of the scalar dissipation rate. The results 
obtained from these two methods are as accurate as 
each other and, since the two methods are derived in 
different ways, we can argue that the Wiener filter is 
a credible approach in evaluating the scalar dissipation 
rate in this flow. 


The analytical tractability is an advantage of the 
Wiener filter, as compared to other filtering methods 
that have been used for making scalar dissipation rate 
measurements. Furthermore, it does not rely on ad hoc 
assumptions and is robust in the spectral modelling in¬ 
volved. Also, since it treats explicitly both spatial res¬ 
olution and noise effects, it can point to potential im¬ 
provements in either of these. 

A detailed presentation of the effect of filtering in 
the instantaneous raw images showed that the correc¬ 
tions in the mixture fraction field were less than 1% 
of the measured value for most times ASI and in the 
case of the scalar dissipation rate these were of the or¬ 
der of 40%. The residual noise after the application of 
the filter in the mixture fraction fields was gaussian 
and constant at the smallest scales, which shows that 
the Wiener filter was capable of largely removing small 
scale measurement noise rather than flow information. 


The procedure outlined in section 4.1 is, also, a way 
to approach the application of the Wiener filter in other 
flow configurations, e.g. flows with higher Reynolds num¬ 


bers or flows with combustion. In all situations the re¬ 
sults of the filtering operation should be tested and 
verified against the filtering assumptions and, if pos¬ 
sible, independently of the applied filter, similarly to 
the present comparison with the method of |Tanaka and| 


Eaton (2007). 


High values of the scalar dissipation rate were shown 
to be concentrated at the jet boundary and at random 
regions within the jet body. Before filtering, noise domi¬ 
nated the instantaneous scalar dissipation distribution, 
at least qualitatively, which was removed after the ap¬ 
plication of the Wiener filter. The instantaneous gradi¬ 
ents were shown to remain after filtering and regions of 
low dissipation values were found, as should be expected 
of the instantaneous scalar dissipation distributions. 

Finally, use of a model spectrum for turbulent flows 
with high Taylor-scale Reynolds number and the present 
experimental conditions, revealed that the Wiener fil¬ 
ter produces errors similar to the errors found in the 
measurements. The separation between the Batchelor 
length scale and fall-off of the PSF and, also, the signal- 
to-noise ratio are the most important parameters in de¬ 
termining the effectiveness of the Wiener filter. 
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